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On Fractional Tempered Stable Motion 


C. Houdre* and R. Kawai^ 
February 28, 2005 


Abstract 

Fractional tempered stable motion (fTSm) is defined and studied. FTSm has 
the same covariance structure as fractional Brownian motion, while having tails 
heavier than Gaussian but lighter than stable. Moreover, in short time it is close to 
fractional stable Levy motion, while it is approximately fractional Brownian motion 
in long time. A series representation of fTSm is derived and used for simulation 
and to study some of its sample path properties. 


1 Introduction 

Fractional Brownian motion (fBm) and its various extensions are not only rich mathe¬ 
matical objects but have also been extensively used in application to model asset price 
dynamics, data traffic in telecommunication network, daily hydrological series, and turbu¬ 
lence, to mention but a few topics. We recall that standard fBm {B^ : t G M}, H G (0,1] 
is a centered Gaussian process with continuous paths and with the following covariance 
structure; 

Cov(i?f, 1 + \s\^^ -\t- , t,seR. (1) 

However, modeling drawbacks of fBm and of some of its extensions have also been dis¬ 
cussed in the literature. For example, although Gaussianity provides analytical tractabil- 
ity, its light tails are often inadequate for modeling higher variability observed in various 
natural phenomena. On the other hand, non-Gaussian stable generalizations immedi¬ 
ately lead to infinite second moment and to a lack of closed form for the density, resulting 
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in significant analytical difficulties. Moreover, the selfsimilar and stationary increments 
properties of fBm are sometimes unrealistic in practical modeling. 

In order to remove these drawbacks, we introduce and study fractional tempered stable 
motion (fTSm), which has the following properties: 

(i) Its marginals have tails heavier than Gaussian but lighter than (non-Gaussian) 
stable (Proposition 1231). 

(ii) It has the same second order structure as fBm fProposition Id.lj) . 

(iii) In long time it behaves like fBm, while in short time it is more akin to fractional 
stable motions, fTheorem Ib.dj) . 

We present its series representation fProposition l4.1|l . which is potentially useful for simu¬ 
lation and which is also used to study sample path properties. When He (l/a, l/a-|-l/2), 
it has a.s. Holder continuous sample paths with exponent (0, if — l/a) fProposition 15.2j) 
and is not a semimartingale (Proposition 12121) • contrast, the sample paths of fTSm 
become nowhere bounded as soon as H < l/a fProposition 15.1|1 . 

Let us close this section by introducing some notations and dehnitions which will be 
used throughout the text. is the f-dimensional Euclidean space with the norm || ■ ||. 
Mg := M'^\{0}, and i3(Mg) is the Borel cr-£eld of Mg. (G, JF, P) is our underlying probability 

space. fl(H) is the law of the random vector Y, while “=” denotes equality in law, or 
equality of the hnite dimensional distributions when stochastic processes are considered. 
Similarly, is used for convergence in law, or of the hnite dimensional distributions, 
while denotes the weak convergence of stochastic processes in the space £*([0, cx)), M) 
of cadlag functions from [0, cx)) into M equipped with the Skorohod topology. G([0, cx)), M) 
is the space of continuous functions from [0, cxd) to M endowed with the uniform metric. 
A sequence of stochastic processes {Xf : t > 0}„gN in C'([0, cxd),M) is said to be tight if 
for each compact set K E [0, oo) and each e > 0, 

limlimsup P [ sup |X” - > e ) =0. 

^^0 n—>oo J 

A sequence of stochastic processes {Xf : t > is said to converge uniformly on 

compacts in probability (ucp) to a stochastic process {Xt : t > 0} if for each compact set 
K E [0, CX)) and each e > 0, 

lim P (sup \X]f - Xt\ > e 1 = 0. 
n^co J 

This last convergence will be denoted by “X" X”. Finally, In’*'a = lira if a > 1. 

In"'' a = 0 otherwise. 

2 Definition of Fractional Tempered Stable Motion 

We begin by briehy reviewing tempered stable distributions and processes (Rosihski [13j i 
which are building blocks for the fTSm dehned below. Let /i be an inhnitely divisible 
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probability measure, without Gaussian component, on Then, /x is called tempered 
stable if its Levy measure has the form 

p poo 

z/(5) = / / lB{sx)s~'^~^e~'^dsp{dx), B E i3(Ro), 

JRg Jo 

where a E (0, 2) and where p, the inner measure^ is such that 



x\pp{dx) < oo. 


( 2 ) 


The two parameters a and p above uniquely identify the Levy measure of tempered stable 
distributions. Under the additional condition 


/ fud W^WpiJ^) < oo, if a e (0,1), 

\/jRd ||a;||(l + \\x\\)p{dx) < oo, if a = 1, 

the characteristic function of p has a closed form expression given by 


p{y) = exp 


for some 6 G R'^ and where, for s E 


I (l>a{{y,x))p{dx) 


s = 


r(— a)((l — is)" — 1 + ias), if a G (0, 1) U (1, 2), 
(1 — is) ln(l — is) + is, if a = 1. 


(3) 


(4) 


Below, we write p ~ TS{a,p;b) if p is given by 0 and denote by {X™ : t > 0} a 
tempered stable Levy process in R such that ~ TS{a, p-,b). Setting 6 = 0 gives 
= 0 for every f > 0, and then {Xj^ : t > 0} is a martingale. 

From now on, we always assume that for any t > 0, 


E[XP] = 0, 


and further that 


\x\‘^p{dx) < oo. 


(5) 


so that for any t > 0, E[(X^^'^)^] < +cx). 

We will dehne fTSm as a process of stochastic integral with respect to the tempered 
stable process, i.e., f(t, s)dX'f^ : t E T}, where f : T x5—s>Risa deterministic 
function. Various such representations of fBm has been introduced in the literature. 
The moving-average representation (Mandelbrot and Van Ness j7]) and the harmonizable 
representation are the most commonly used. These have also been extended to non- 
Gaussian stable marginals. (See Ghapter 7 of Samorodnitsky and Taqqu [IH].) A lesser 
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known representation of fBm involves a Volterra kernel and is due to Decreusefond and 
Ustiinel jS|. Recall that a Volterra kernel is a function R': MxM —>■ M such that K{t, s) = 0 
for s > t. In the present paper, we will use the Volterra kernel KH,a '■ [0, cxo) x [0, cx)) —*• 
[0,oo), given by 




^H,o. 






^1/a-H 





( 6 ) 


where H G (1/a 


1/2,1/a + 1/2), a G (0, 2), and 


/G(l -2G)r(l/2-G)y/^ 
Vr( 2 - 2 G)r(G+ 1 / 2 )) ’ 


with (throughout) G := H — 1/a + 1/2. Clearly, Ki/o,,a{t,s) = l[o, 4 ](s). Note that when 
H G (1/a, 1/a + 1/2), we also have 

s) = chAH - l/a)s^/"“-^ [ A - l[o.4](s). (7) 

J S 

Despite its complex structure, there are two main advantages to using this kernel as an 
integrand: 


(i) It is dehned only on [0,f], f > 0, while the domain of dehnition of the moving- 

average kernel {t — s)+ — (—s)+ is the whole real line. (A moving-average 

fractional Levy motion has recently been studied in Benassi et al. ^J.) In gen¬ 
eral, however, it is impossible to generate background driving stochastic processes 
dehned on M. 

(ii) It can be written as a Riemann-Louiville fractional integral, whose inverse function 
has a closed form which is also a Riemann-Louiville fractional derivative. This is 
important for the prediction of the sample paths of fTSm, problem of consequence 
in hnancial modeling. This will be presented in a subsequent paper |3]. 


Below, we derive a necessary and sufficient condition on H so that for each f > 0 and 
p>2, the kernel is in L^([0,t]). 

Lemma 2.1. Let t > 0, let a E (0,2), and let p > 2. Then, KhA^e) ^ -^^([0,^]) if and 
only if H G (1/a — 1/p, 1/a-|-1/p). In particular, KhA'^e) ^-^^([0,t]). Moreover, when 

KnA'^y ■) ^ -^P([0, f\). 


where 

^H,a,p 



"S)^cis p 


^p{H—l/a)+l 


( 8 ) 




w 




{w — vA 


dv. 


(9) 
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Proof. The case H = 1/a is trivial since Ku^aifis) = l[o,t](s). If if > 1/a, then 
KH,a{'t,s) > 0 , Kh,0(^,3) is decreasing in s, and K^,0(^,3) (j'gi/a ^ as s I 0 for 
some constant C'. Hence, KH,a{t,-) G LP{[ 0 ,t]) if and only if p(l/a — H) > —1, i.e., 
H < 1/a + 1/p. When H < 1/a, KH,a{.-,3) explodes at s = 0 and s = t. In fact, 
Kh,0(^,3) ~ as s J. 0 and KH,a{t, 3 ) ~ C"'{t — as s t t for some con¬ 

stants C" and C". Thus, Ku^aii,-) ^ -^^([ 0 ;^]) if and only if p{H — 1/a) > —1, i.e., 
H > 1/a — 1/p, which proves the hrst claim. The last claim follows from elementary 
computation. □ 

Remark 2.2. Above, we only considered the case p >2 because of the moment condition 
(0 and since H G (1/a — l/2,l/a-|-l/2) in (jHl). However, it is easily seen that the results 
above remain true for arbitrary p > 0, provided that the kernel is defined for arbitrary 
H. In particular, we have Kuwait, ■) G T"([0, t]) since (0, 2/a) D (1/a — 1/2,1/a -|- 1/2). 

Let us state two other known properties of the kernel. The proof of (ii) below can be 
found in, e.g., Decreusefond and Ustiinel jH], or Nualart jS], while (i) is immediate. 

Lemma 2.3. (i) For each h > 0, 


KH,a{ht,s) = ^/°‘KH,a{t,3/h). 

(ii) For t,s > 0, 



where G = H — 1/a + 1/2. 


We are now in a position to dehne fTSm. 

Definition 2.4. Fractional tempered stable motion {Lf : t > 0} m M is given by 



( 10 ) 


where the integral is defined in the L^(r2, JF, P) -sense. 

The integral above is well dehned by the moment condition ®, Lemma 12.31 (ii) and 
the help of the Wiener-Ito isometry. (See also the proof of Proposition 13.II below.! For 
convenience, we will henceforth write t > 0} ~ fTSm{H,a, p) when {L(f : t > 0} 

is dehned as m- We note that = X'^^, which is a Levy process, since Ki/a,a{t, 3 ) = 
l[0,t]('S)- 

The following is an important result on the marginals of fTSm. 

Proposition 2.5. The finite dimensional distributions of fTSm are tempered stable with 
finite second moment. 
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Proof. Let A; G M. It suffices to show that for any real sequence and any non¬ 
negative nondecreasing real sequence the random variable Yli=i is tempered 

stable. First, observe that 


^ /‘ft, / k 

(e- 

i=i -^0 


dxj^. 

Then, by Proposition 35 of Rocha-Arteaga and Sato mil, we get 


E[( 





k 

aiKH,a{ths))p{dx)ds 

1=1 


exp 


(j)a{,yx)r]{dx) , 


where (fa is given by 0 and where p = M o J with M{dx, ds) = p{dx)ds and 


J{B) = I (x,s) e Mo X [0,4] : x'^aiKn^aiU, s) G R > , Be i3(Mo). 


i=l 


The measure p is well dehned as an inner measure with hnite second moment since for 
each i, Ku^aiU:') ^ -^^([04*]) 


\x\‘^p{dx) = / \x\‘^p{dx) 


J ^ Q,(t^, 


ds < 00 , 


which concludes the proof. □ 

It is worth noting the one dimensional marginal result as a corollary. 

Corollary 2.6. Let {L^ : t > 0} ~ fTSm{H,a, p) and let (fa be given by (H- For each 
t > 0, 


] = exp 



(fa{yxKH,a(t, s))p{dx)ds 


IJo 


( 11 ) 


and thus 

Lf ~T5(a,r/p0) (12) 

where pt = M o with M{dx,ds) = p{dx)ds and Jt{B) = {(a:,s) G Mq x [0,t] : 
xKH,a{t,s) G B}, B G -B(Mo). 


3 Covariance Structure and Long-range Dependence 

Let us hrst describe the covariance structure of fTSm. 
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Proposition 3.1. Let {Lj : t > 0} ~ fTSm{H,a, p). Then, 'E,[Lf] = 0, and 


Cov(Lf, Lf) = ^ - It - mx^n s,t>0. (13) 

Proof. Recall that {Xf^ : t > 0} is a square-integrable centered martingale. The hrst 
claim is thus trivial. For the second claim, observe that for s G [0,t], 

Cov(Lf,Lf) = E[LfLf] 


= E 

= E 


KH,a{t,u)dXl^ I KH,a{s,u)dX, 


TS 

u 


Uo 


rsAt 


fsAt 


KH,ait,u)dX.^ 


TS 


KH,a{s,u)dX,^ 


TS 


L-'O 


rsAt 


= E[(X/'^)^] / KH,a{t,u)KH,a{s,u)du, 

Jo 

where the last equality holds by the Wiener-Ito isometry. Lemma 12.31 (ii) then gives the 
result. □ 

Let us state some immediate consequences of the previous result. 

Corollary 3.2. Let {L^ : t > 0} ~ fTSm{H,a, p). For each t > 0 and each h > 0, 

= (14) 

and for s,t > 0, 


- L^Y] = = \t - srE[{xrY]. 


H n 2 i 


2Gti 


-rs\2i 


(15) 


The property (d is sometimes called second-order self similarity. Moreover, m says 
that fTSm has second-order stationary increments, which clearly implies its continuity in 
probability. 

We are now in a position to discuss the long-range dependence of fTSm. The dehnition 
of long-range dependence is often ambiguous and varies among authors. In the present 
paper, we will follow Samorodnitsky and Taqqu d; the increments of a second-order 
stochastic process {Xt : t > 0} exhibit long-range dependence if for each h > 0, 

OO 

^ ^ I Cov(^/2, ^nh ^(n—l)/i) | OO^ 

n=l 

or short-range dependence^ if each /z > 0, 

OO 

I Cov(X/^. — Xq, Xnh — X(^ri-l)h) \ < OO. 

n=l 
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Proposition 3.3. The increments of fTSm exhibit long-range dependence when H G 
(l/o, l/a+1/2), and short-range dependence when H G (l/o — 1/2, l/o]. 

Proof. By Lemma (3.11 we have for each h > 0, 

Cov(L^, - Lf) = + h/tf^ - 2 + (1 - h/tf^) 

_ - l)h2, 

as t —*• oo. The claim then holds since 2{G — 1) > —1 for H G (l/o, l/o + 1/2), while 
2(G-1) <-1 for i/G (l/o-1/2,1/a]. □ 

In relation to the second moment, we will consider higher moments of fTSm. 

Proposition 3.4. Let {Lf : t > 0} ~ fTSm{H,a, p). Then, for each p > 2 and each 
t > 0, E[|L/^|p] < oo if and only if H E (1/a — 1/p, 1/a + 1/p) and \x\^p{dx) < oo. 

Proof. By Corollary 12.14 for each t > 0, C{L/f) is tempered stable. By Proposition 2.3 
(hi) of Rosihski d, E[|Lf |P] < cx) if and only if \x\Ppt{dx) < oo, where pt is the 
inner measure of L/f given as in Corollary 12.hi that is. 


{\x\KH^ait, s))^p{dx)ds < oo. 
The left hand side of the above can be decomposed into two terms; 






<|x|< 




VI 


xf p{dx)ds. 


and ^ 

f KHa{t,s)'^ f \x\^p{dx)ds. 

The hrst term is equivalent to Kuwait, s)'^ds due to the moment condition (jS)) on p, while 
the second terms is clearly equivalent to KH,a{t, s)'^ds \x\^p{dx). Then, LemmaEUl 
concludes the proof. □ 

Remark 3.5. It is shown in Proposition 2.3 (iv) of Rosihski Hg that a tempered stable 
distribution has exponential moment of certain order if and only if its inner measure has 
a compact support. Unfortunately, the tempered stable marginal of fTSm cannot have 
exponential moment since, with the notation of the preceding proposition, for any e > 0, 


Pt{{x G Mo : \x\ > e}) = 


'\x\KH,a{t,s)>e 


p{dx)ds > 0, 


due to the unboundedness of the kernel KH,a- 





4 Series Representation 


In this section, we derive a series representation of fTSm, which is inherited from the one 
of tempered stable processes obtained by Rosihski jHj. This representation can also be 
used for simulation. (See Figure HI) Moreover, we will make use of its structure to derive 
some sample path properties in Sectional 

Let {rj}j>i be arrival times of a standard Poisson process, let {Ei\i^i be a sequence 
of iid exponential random variables with parameter 1 , let {Ui\i>i be a sequence of iid 
uniform random variables on [ 0 , 1 ], let be a sequence of iid random variables in 

Mq with common distribution 

\x\°'p{dx) 

miypY 


and let be a sequence of iid uniform random variables on [ 0 ,T]. Also let m{p)'^, 

k', and zt be constants given by m{p)^ = |x|“p(da;), k' = m{p)~'^ ^^^x\x\°'~^p{dx), 

and 


m{p){a/T) ^ + |r(l - a)| if a 7 ^ 1, 

(In(m(p)T) + 27 ) f^^xp(dx) - f^^xln \x\p{dx), if a = 1, 


where C, denotes the Riemann zeta function and 7 (= 0.5772...) is the Euler constant. 
Then, Theorem 5.4 of Rosihski HI] tells us that 


E 




— IjOL 




Ye 

1^*1 


l(Ti < t) 


m{p) 



+ zrt 


converges a.s. uniformly in t G [0,T] to a tempered stable process with TS'(a,p;0). This 
series representation can easily be extended to fTSm as follows. 

Proposition 4 . 1 . Let {L^ : f > 0} fTSm{H,a, p) and let T >0, Then, 


{if :«e|0,r]} 


E 

2 = 1 


m{p) 


aTi 


— llOL 


^E,uYm 1 


Vi 


Mp) ( Y 


T 

-l/« fH-l/a+l 




+ ZtCh, a ,: t E [0,T] ^(lO) 


where CH,a,i is the constant defined by If p is symmetric, then 


{if :«e|o,r]}^l E(™W 


2=1 


oTj 

~¥ 


— Ija 


LE,uE\V,\ I Y«-aXTij : ( € [0,T] 7 


V, 


Moreover, if H E [1/a, 1/a + 1/2), then the series converges almost surely uniformly in t 
to fTSm{H, a, p). 
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Proof. We will only consider the asymmetric case. By arguments as in Theorem 5.4 of 
Rosihski El, we get 


E 

2 = 1 


C(i\ 4.H-l/a+l 

m{p) ( — j k'CH,a,i^ - c,{T)E[KHAt,T,)] 


— ZTCH,a,lt^ 


uniformly in t, where 

AT) := / E 


'i-l 


m{p) 


ar 


— Ija 


AE,Ul^‘^\V,\] ^ 


l^il 


dr. 


Hence, the right hand side of (ITBl) can be rewritten 


as 


E 

2=1 


m{p) 


aTi 

T 


— 1/a 


AE,UI^^\V,\ ) ^^KH,aAT^) - AT)E[KHAtAi)] 


Vi 


:= Zi 


(17) 

Next, we need to analyze the hnite dimensional distributions of {Zt : t G [0,T]}. Let 
k eN, let {aj}Ai be a real sequence, and let {tj}Ai be a nondecreasing sequence taking 
values in [0,T]. We will show that the random variable bas the same law as 

■ III view of Proposition 12.5L we have 


E 


= exp 

pT p ^ 

/ / (fajyx Y. ajEnAtj, s))p{dx)ds 




Jo J Mq j=l 


where (fa is given by 0- Also, observe that 


3=^ 


2 = 1 


m{p) 


aTi 

~Y 


— 1/a 


A 


EiUAm ) 


V. 


-Ci{T)E 


k 

^ ^ [tj ; El 

i=i 


This series representation is induced by the Levy measure 
u{B) 



oo nl noc 



Is(i7(r/T, M, s, x) Y/j o>jKH,Aji v))drdue ‘^dspi{dx) — 

pi p poo ^ 

/ / / 1b{sx E '^dsp{dx)dv, 

Jo Jmo Jo -^1 


^ 0 J Mo ^0 ^0 ^0 

rT 


where H(r,u, s,x) = {m{p){ar) A s-u^/"|a;|)a;/|a;| and pi{dx) = m(p) "|x|“p(da;). In 
fact, the measure v is well dehned as a Levy measure since KH,a{tj^') ^ E^([0,tj]) for 
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each j. Therefore, by Theorem 4.1 (B) of Rosihski |i;fj . 




exp 

f _ 1 _ iyz)v{dz) 

.J Rq 


exp 

' fT r ^ 

/ / (t)a{yx^ajKH,o 

■JO -'Ro 

ftj, s))p{dx)ds 


which proves the equality of all hnite dimensional distributions. 

Next, let H E [l/a, l/a + 1/2). Dehne, for t G [0, T] and s G [0, cxd). 


Zt,s ■— 


E 

{i:ri<s} 


m{p) 


aVi 

T 


-1/ol 


AE,a 


1 /oL 


K\ I 


-Ci(r)E|A-„,„(«,ri)] 


Notice that : t G [0,T], s G [0, cx))} has independent increments in s {not in t, of 
course.) Arguments as in the proof of Theorem 5.1 of jT^] give the a.s. convergence of 
the series uniformly on [0,T]. □ 


5 Sample Path Properties 

In this section, we investigate sample path properties of fTSm. Let us begin with the case 
H G (1/a - 1/2,1/a). 

Proposition 5.1. When H G (1/a —1/2,1/a), fTSm is a.s. unbounded on every interval 
of positive length. 

Proof. Let T > 0. For each t G (0,T], lirngp KH,a{t, s) = +cxo. Hence, in the series 
representation given in Proposition 14.IL KH,a{Ti,Ti) = +oo, for all i G N, and so none of 
the summands are well dehned. This implies that supjg[o,T] ^ 

Unfortunately, the above sample path property makes fTSm with short-range de¬ 
pendence of little practical use. We notice that for any H G (1/a — 1/2,1/a -f 1/2), 
EH,a{t,0) = +00. But, this turns out to be irrelevant to the unboundedness of sample 
paths since Tj 7^ 0 a.s., i G N. 

FTSm with long-range dependence has better sample path properties. In particular, 
it has a Holder continuous version with exponent 7 G (0, if — 1/a). 

Proposition 5.2. If H E (1/a, 1/a-f 1/2), there exists a continuous modification of 
fTSm, which is a.s. locally Holder continuous with exponent 7 for every 7 G (0, if — 1/a). 

Proof. By Corollary 1,4.21 we have ^[\Lfi — = \t — sp‘^E[(X/’'^)^]. If if > 1/a, then 

2G > 1, and thus the Kolmogorov-Centsov Theorem (see, for example. Theorem 3.23 of 
Kallenberg j3]) directly applies, giving the result. □ 
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We will henceforth always assume that when H G (1/a, l/a + 1/2), we are using such 
a continuous version of fTSm. 

Proposition 5.3. Let {L^ : t > 0} ~ fTSm{H, a, p) with H G (1/a, 1/a + 1/2). 

(^) 

fiV-l 1 


(a) FTSm is a.s. 
(Hi) FTSm is not 

Proof, (i) Immediate 


lim E 

N^OO 



= 0 . 


_n=0 

of infinite variation on every interval of positive length, 
semimartingale. 

from 


E 


TAf-l 


E \tH - tH 

N 

,n=0 


2 




(ii) Let T > 0. For each s G [0,T], we have 

\KH,aitl,s) - KH,ait2,s)\ 


lim sup 


1^1 ~ ^2 


+ 00 , 


which implies that for each s G [0,T], KH,a{-,s) is of inhnite variation. By Theorem 4 
of Rosihski m and a symmetrization argument given there, fTSm is of inhnite variation 
with positive probability. FTSm is selfdecomposable and hence by Corollary 3 of Rosihski 
[12], it obeys a zero-one law. This gives the result. 

(hi) The convergence in (i) implies convergence in probability, and together with (ii), 
the claim follows from the same argument in Lin [Uj. □ 

Remark 5.4. In view of (iii) above, the stochastic integration for fTSm cannot be dehned 
in the classical semimartingale sense. However, a slight modihcation of the kernel 
induces a corresponding semimartingale, which can be arbitrarily close to fTSm. Recall 
the simpler expression of the kernel for H G (1/a, l/a-t-1/2) given by (|7|). Observe that 


s) = CH,a{H - l/a){t - s)^ ^ l[0,p(s), 


and that 


Therefore, 


/■* d 

KH,a{t,s)= / —KH,a{u,s)dul[o^t){s). 


if = / KhAV s)dxp = / (/ «)*) 


(18) 


If the two integrals could be interchanged, then fTSm would be of hnite variation, i.e., it 
would be a semimartingale; we have seen that this is not the case. On the other hand. 
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the integrability condition of the stochastic Fubini’s theorem (Theorem 46 of Protter P) 
can be achieved by slightly modifying the kernel Set 





H-l/a-l 


and so we have 


n e M, 
(19) 


Wt 




chAH - 1/a) 



\ H—I/ol—1 

) 



l[o,t](s)- 


The integrability condition is then satished; for every u G [0,f], 


d_ 

du 


= clAH - 1/af I 
< c]jAH - l/oifil - - l/a))-^n-2(-f^-V“-i)^i < oo, 


^ \2{H-l/a-l) ,^.2(H-l/a) 


M H- S 

n 


ds 


and thus the stochastic Fubini’s theorem applies. Therefore, we get 

J‘ Kh Jt, s)dxr = j‘ du, 

which is clearly of hnite variation. For hnancial modeling, it is of interest to further 
modify the above to get an inhnite variation semimartingale. This can be done as follows. 
For e > 0, set s) := s) + e. Since = §1 ^hA^^ ■®)’ stochastic 

Fubini’s theorem again applies and thus we get 

[‘K"„At,s)dxr = [\c+Kijt,s))dxr 

Jo Jo 

which exactly follows the dehnition of the canonical decomposition of semimartingales, 
i.e. a martingale plus a hnite variation process. 


6 Short and Long Time Behavior 

In this section, we obtain the short and long time behavior of fTSm, which are also 
inherited from the background driving tempered stable processes obtained by Rosihski 
jl4j . In short time, fTSm is asymptotically fractional stable motion (to be dehned below), 
while in long time it is approximately fBm. 
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Let us begin by briefly reviewing the corresponding behaviors of tempered stable 
processes, proved in Theorem 3.1 of m. 


(i) Short time behavior . Let {Xf^ : t > 0} ~ TS'(a, p; 0) and let 

{ hT{l-a) j^^xp{dx), i/ae(0,1), 

-(1 + Inh)i/a = 1 , 

0 , */ae(l, 2). 

Then, 

: t > 0} 4 {X“ : t > 0}, as h ^ 0, 

where {X“ : t > 0} is an a-stable process in M such that 




= exp 


t / pa{yx)p{dx) 


where 


( 20 ) 


Pa{s) = 


—r(—a) cos ^|s|“(l — itan ^sgn(s)), if a G (0,1) U (1, 2), 


— (||s| + is In |s|) + is, 


if a = 1. 


( 21 ) 


(ii) Long time behavior. Let {Xj^ : t > 0} ~ TS{a, p] 0). Then, 

: t > 0} 4 {cWt : t > 0}, as h ^ oo, 
where {Wt : t > 0} is a standard (centered) Brownian motion and 

= r(2 — a) [ x^p{dx). (22) 

J Mq 

We will say that the limiting stable process {X" : t > 0} given above is associated to 
the tempered stable process : t > 0}. 

Let us now define fractional stable motions (fSm), which turns out to be a short time 
limit of fTSm. Below, for a 7^ 1, the integral in (I23|l is well defined in probability since 
Xn^aifi') £ -^“([Oj^]); t > 0. On the other hand, when a = 1, the extra symmetry 
assumption on p also ensures that it is well defined in probability. (See Remark 12.21 and 
Samorodnitsky and Taqqu uni-) 

Definition 6.1. Let {L^ : t >0} fTSm{H,a, p), where when a = 1, p is additionally 
assumed to be symmetric. Fractional stable motion (fSm) {L^’°‘ : t > 0} associated to 
fTSm {L^ : t > 0} is given via 

L^’^-.= [ KHAt,s)dX:, t>0. (23) 

Jo 
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Let us derive some basic properties of fSm. 

Lemma 6.2. Let t > 0} be fSm associated to {Lf : t > 0} ~ fTSm{H, a, p). 

(i) The finite dimensional distributions of fSm are stable. 

(a) For each t > 0, the characteristic function of is given by 


E[e 


T ti ,a 

— 


= exp 


C 


H^a^a 


aH 


'Pc,(yx)p{dx) 


( 24 ) 


where CH,a,a is the constant given by and where 


(Pa(s) = 


(Pa(s), i/a e (0,1) U (1,2), 


2 


if a = 1, 


H.a 

s • 


where ipa is defined by (HI. 

(Hi) FSm is selfsimilar; {h~^L^f : t > 0} = t > 0}. 

(iv) FSm has (strictly) stationary increments; for each t > s > 0, = LfS 

(v) When H G (1/a, 1/a + 1/2), fSm has a continuous version, while when H G 
( 1 /a — 1 / 2 , 1 /a), it is unbounded on every interval of positive length. 

(vi) When H = 1/a, fSm is an a-stable (Levy) process. 

(vii) With the notation in Section 4, if ct G (1,2), then 


{Lf’“:tG [0,T]} 


E 

2=1 


m{p) 


aTi 

T 


-V« y, 

/ a \ —1/“ 

+m(p) [fj ((l/a)k'CH,a, 


1/a+l 



while if a & ( 0 , 1 ], or if a E ( 1 , 2 ) and p is symmetric, then 


{L«'’‘:t€lO,T]}£ {J2 Mp) 


2=1 


aTi 


— Ifa 


Yf 

1^*1 


Kuwait, Ti) -. t E [0, T] 


Moreover, if LI E [1/a, 1/a + 1/2), a G (0,2), the above series converges almost surely 
uniformly in t E [0, T] to a version of {L(^’°‘ : t E [0, T]}. 
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Proof. With the notation of Proposition 12.51 for a ^ 1, 


E[e' 


E k T H ,OL 


= exp 


— r(— a) cos 


TTP 




^0 •/Rq 


y^^aiKH,a{ti, s)yx 


Z =1 


/, . Tra ( 

1 1 — ■? tan —sgn 1 

k 

“ P 1 

yx I I p{dx)ds 

\ ^ \ 

2=1 

) ) 


= exp 


■p/ \ 

f (—a) cos — 




ds 


i=l 

/ \yx\°‘(l — ita.n—sgn{yx)j p{dx) 
Jro ^ 2 / 

When a = 1, the symmetry of p yields = 0, and so 


E[e 


• V~' ^ T 1 

wEi^iCiiP. j ^ 


= exp 



TT 

2 


'^aiKH,i{ti,s)yx 


2 = 1 


p{dx)ds 


TT 

2 




2=1 


ds / \yx\p{dx) 


For each a e (0,2), XiLi •) ^ L"([0,4]) since aiKn^aitir) ^ -^“([0,tj]) for each 

i. Thus, (i) holds. Clearly, (ii) is a direct consequence of (i). (iii) and (iv) follow from 
the selfsimilarity and stationary increments properties of : t > 0 } with the help 

of Lemma ESI By (iii) and (iv), = \t — s|^^F'[|Lf’“|^], 0 < p < a. 

Hence, when H G (1/a, 1/a + 1/2), the continuity follows from the Kolmogorov-Centsov 
Theorem. Next, let H G (1/a — 1/2,1/a). For each T > 0, supjg[o,T] l-^ir,a(Cs)| = +cxd, 
s G [0,T]. The nowhere boundedness follows from (i) and Theorem 4 of Rosihski ^T] as 
well as the zero-one law for stable processes and a symmetrization argument given there. 
Hence, (v) holds, (vi) is immediate from Ki/a^ait,s) = l[o,t](s). Finally, (vii) follows from 
arguments as in Proposition 5.5 of Rosihski ^3] and from Proposition 14.11 □ 

We will henceforth always assume that when H G (1/a, 1/a-|- 1/2), we are using a 
continuous version of fSm. 


Remark 6.3. Many extensions of fBm are available in the stable literature; for example, 
linear fractional stable motion, log-fractional stable motion, harmonizable fractional stable 
motion. (See, e.g., Samorodnitsky and Taqqu [lEl-) These various extensions are not 
necessarily identical in law since their marginals are determined by kernels in the stochastic 
integral representations. Indeed, fSm dehned above is still different from any of them. 


We are now in a position to present the main result of this section. 
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Theorem 6.4. Let {L^ : t > 0} ~ fTSm{H,a, p) with H ^ 1/a. 
(i) Short time behavior : Let 


b 


r(i*/ae(o, 1), 

0, tfae[l,2), 


and let kt = Kuwait■, s)ds. Then, 


{h ^L^^ + h^ : t > 0} ^ > 0} as h ^ 0, 


where {L^’'^ : t > 0} is fSm associated to {Lf : t > 0}. 


(ii) Long time behavior: 


{h t > 0 } {cBf : t > 0 } as h ^ oo, 

where {B^ : t > 0} is a standard fBm with parameter G and where c is the constant given 

by (d^. 

(Hi) When H G (1/a, 1/a + 1/2), the convergence in (i) and (ii) can he strengthened 
to the weak convergence in (^([O, cx)), M). 

Proof of (i) and (ii). (i) Observe that for each t > 0, 

h-^L((t + h^-^l^bkt = [ KnAt, s)h-^/^d{X^f + bhs) := 

Jo 

It thus suffices to show that for any real sequence {a^^i and nonnegative nondecreas¬ 
ing real sequence k E N, the random variable converges in law to 

Y)’i=i as h ^ 0. Since 


EW = I (E' 

i=l >^0 \i=l 


iKH,aAw) I h + bhs), 


we have by Proposition 12.51 that 



= exp 

rtk r ^ ^ 

/ / aiKnAti: s))p{dx)ds 

*■ 


Jo JMq i=\ 


where 

( r(—a)((l — is)“ — 1), if a e (0,1), 

V’o('S) = I ln(l -|- s^) — s tan“^ s, if a = 1, 

[r(—a)((l — is)°‘ — 1 -|- ias), if a G (1, 2). 


(25) 


(26) 


Note that "i/i is obtained via the symmetry of p. (See Proposition 2.8 of Rosihski ^1].) We 
then want to show that (ESI) tends to the characteristic function of the random variable 
as h ^ 0 . 
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The proof of Theorem 3.1 (i) of Rosihski HU shows that for a 7 ^ 1 , 

lim h'ipa{h~^^°'s) = 
h^O 


where (pa is given by m and 

< Za\s\°‘, 

where Za is some constant depending only on a{^ 1). When a = 1, we have 

TT 

lim hil)i{h~^s) =-|s|, 

ft —>0 2 

and the uniform boundedness (in h > 0 ) of \h'ipi{h~^s)\ can be shown as 

|h' 0 i(h“^s)| < |hln Vl + + |stan“^(h“^s)| 

< |hln(l + h“^|s|)| + ^|s| 



Clearly, Yl\=i (^iKH,a{U-, s) G T"([0, tk]) since aiKH,a{U-, s) are in T"([0, U]). Together with 
the moment condition m on p, the passage to the limit in (123) is justihed. Hence, 


lim E 

ft^O 




rtk j- 


= exp 

/ / ipaiyx} ^aiKH,aiti,s))p{dx)ds 

*■ 


Jo J Mo i=\ 


which is the characteristic function of X]i=i 
(ii) We have that for each h > 0, 


Cov(h-«L^„ h-^Lg) = ^ - (t - EiiXrn s e [ 0 ,t]. 


Hence, for the convergence of all hnite dimensional distributions, we only need to show 
that the marginal law at any time of {h~^Lg : f > 0} converges to Gaussian. Without 
loss of generality, let t = 1. By Lemma f2.31 (i). 


E[e*^^ = gxp 

/ / h„(h ^yxKH,a{h,s))p{dx)ds 

_J 0 J Mq 

= exp 

[ [ 'da{h~^^'^yxKH,a{^,s))p{dx)ds 

_J 0 J Mq 


where 'da{u) = — 1 — ius)s “ ^ds. As in the proof of Theorem 3.1 (ii) of 

Rosihski HI], h follows that 


\do,{yxKH,a{^,s))\ < {yxfT{2 


a) 



KH,a{^,sYds, 
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which justihes the passage to the limit below 


lim 

h^oo 


exp 


/ KHAhsfdsr{2 

^ 7o 



This shows the convergence to a Gaussian law, which concludes the proof. □ 

To proof (hi), let us present a technical lemma. 

Lemma 6.5. Let {X" : t > cind : t > 0}„gN sequences of stochastic 

processes in G([0, cx)),M). If the sequence {Xf : t > is tight and if for each n eN, 

there exists a version Y^ofY^ defined on the same probability space as X'^ and such that 
the sequence {Yf' — Xf : t > 0} converges ucp to zero, then the sequence {Yf' : t > OjneN 
is tight. 


Proof. For each compact set K G [0, cx)) and for each 5 > 0, we have 


E 


sup 

t^s^K,\t—s\<5 








< E 

sup 

t,sGK,\t—s\<6 

yn _ xn 

A 1 

+ E 

sup |Xf-X;|Al 

t,sGK,\t—s\<S 


+E 


sup 

t,sGK^\t—s\<6 


Xf - YI 


A 1 


2E 

sup 

yn _ xn 

A 1 

+ E 

sup |Xf-X;|Al 


teK 




t,sGK,\t—s\<S 


(27) 


The hrst term in (EH) tends to zero as n — > cx), since Y^ — X"^ 0. The second term also 

tends to zero as n —> cx) and 5 0 by the tightness of {Xf : t > 0}„gN in ^*([0, oo),M). 

The claimed result then follows. □ 


Proof of Theorem (Hi). By, for example. Lemma 16.2, Theorem 16.3 and 16.5 of 
Kallenberg j3], it suffices to show the tightness of the sequences {h ~^+ h^~^/°‘bkt : t > 
0 } (as h I 0) and t > 0} (as h | cx)) in G([0, cxo), R). 

We begin with the short time behavior case. By Lemma 16.21 (iii) and (iv), for each 
P 6 (0,a), 


Ei\h-«L«r - h-XLjirn = (*- orBiiif’-ii, e io.«i. 


(28) 


By Corollary 16.9 of Kallenberg j^, the uniform boundedness in h seen in (I28|l implies 
the tightness of the sequence {h ~^t > 0} in G([0, cx)), R). Hence, by Lemma IH7H1 
it is enough to hnd h~^L^’°‘ and h~^L(f, + h^~^^°‘bk. in G([0,T],R), T > 0, dehned on 
a common probability space and such that h~^L^t°‘ — + h^~^/°'hk.) 0. To 

this end, we make use of their series representations. First, let a G (0,1], or let a G (1, 2) 
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with a symmetric p. With the notation of Section 4, it follows from Proposition 14.11 and 
Lemma o (vii) that the stochastic processes 


h ^^m{p) 

i=l 



-l/a y 

hTi) 


h-^LZ’- 


and 






— 1/a 




A 


KH,a{h-, hTi) 




converge almost surely uniformly on [0,T], respectively, to versions of fSm and fTSm 
h~^L/’°‘ + h^~^^°'bk., dehned on a common probability space by using the common random 
sequences {Ti}i>i, {Vi}i>i, {Ei}i>i and Then, in view of Lemma 0(i), 


h-^L/:^ - + h^-^/^bk.) 


E 

i=\ 


m{p) 


aTi 

T 


— Ifa 


h-^/^EiU] 



(29) 


which clearly converges ucp to zero as h ^ 0. Also, in case a G (1, 2) and p is asymmetric, 
a similar argument yields ()21H1 . Hence, the sequence + h^~^^°'bkt : f > 0} is tight 

in C([0, cx)),M), which concludes the short time behavior case. 

For the long time behavior case. Corollary Id.21 gives for h > 0, 

E[(h-^Lf, - h-^L^)^] = it- s)^^E[{Xr)% (30) 


Again by Corollary 16.9 of Kallenberg j^, the uniform boundedness in h seen in dSHD 
implies the tightness of the sequence : f > 0} in C([0, cx)), M), which completes 

the proof. □ 

Remark 6.6. The short time behavior result can also be seen from the series represen¬ 
tation. For simplicity, consider the symmetric case. With the notations of Theorem 14. IL 
we have by Lemma f2.41 fii 

oo / \ —1/*^ Ty 

^ ^m(p) -^^Kn^ait^Ti) a.s., as h ^ 0. (31) 


Ji 

iv^.l 


KH,a{ht, hTi] 
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We have seen in Lemma o that the marginals of fSm are a-stable and so their 
covariance does not exist. An alternative notion is the one of covariation. For two jointly 
symmetric a-stable random variables X and Y with a > 1; 

r(A,y):=||A||“ + ||y||“-||X-y||“, (32) 

where the norm || ■ || gives the scale of parameter, i.e. for Z ~ S'q,((T, 0, 0), \\Z\\a = cr. More 
generally, one can also dehne the codifference for any jointly inhnitely divisible random 
variables X and Y as 


I{ei,e2-,X,Y) := +lnE[e*®i^] +lnE[e*® 2 V|^ (- 33 ^ 

for 6 * 1 , 6*2 G M. Clearly, is a special case of Let {L^’°‘ : t > 0} be fSm associated 
to {Lf : t > 0} ~ fTSm{H,a, p) where p is symmetric. Then, 

/(l,-liLf’“,Lf'“)=C(f"" + s“"-((-s)"«), S€l0,(], (34) 

for some constant C. In the Gaussian case, the codifference coincides with the covariance. 
For example, for a standard fBm {Bf : t > 0}, 

, Sf) = + 5 ^^ - (f - = Cov(Sf, ), s G [0, f], 

and 

^(BS-1 - Bf, Bf - Bf) = CoMBf, - Bf, Bf - Bf) ~ 

as t —>• cx). See Samorodnitsky and Taqqu M for more details of the covariation and the 
CO difference. 

Interestingly enough, as shown below, the codifference of the increments of fTSm has 
the same order of decay as the covariance (see Proposition 13.3j) . 

Proposition 6. 7. Let {L^ : f > 0} ~ fTSm{H,a, p). Then, 

/(«1. fc; if+i - Lf, Lf - Lf ) ~ G(»,, 

as t —>■ 00 , where 

(F( 6 >i, 6 > 2 ) = r ((1 - ia; 6 ' 2 iFH(l,s))"“^ - 1 ) xs^/""-^dsp(dx). 

F(a)sm( 7 ra) J[o,i]xKo 

Proof. Observe that 

/(0i,02;4^i-Ff,Lf-L^) 

= F(-a) [ ( - (1 - ixiOiiKnAt + 1, s) - s)) + s)))“ 

JrxRq ' 

+ (1 - ixei{KH,ait + 1, s) - KH,a{t, s)))“ 

+ (1 — ix92KH,a{^-i s))°‘ — l^dsp{dx), 
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and that 



as t —i> cx). Hence, for each s > 0, 



+ (1 - + 1, s) - s)))" + (1 - ix92KH,a{^, s))“ - 1 

~ iax9iCH,aS^^^~^ ((1 — ix92KH,a{^-, — l) 


as t —> oo. The resnlt then holds since r(—a) = 


— TT 


□ 


7 Concluding Remarks 

Willinger, Taqqn and Teverovsky ini assert that a nnmerical analysis of stock price time 
series indicates long-range dependence with marginal tails heavier than Ganssian bnt 
lighter than stable. Moreover, it has been known that in shorter time, the asset price 
paths tend to lack higher moments, while they have a Ganssian behavior in long time. 
Indeed, fTSm achieves all those properties. 

In Fignre ^ we give typical sample paths of fTSm and of its backgronnd driving 
tempered stable processes, generated via the series representation presented in Proposition 
o We pnt the inner measnres pi and p 2 as pi{dx) = 5-i,o{dx) -|- 5i,Q{dx) and p 2 {dx) = 
0 .5“"5_o,5(da;)-|-(5i.o((ia;). (Tempered stable Levy processes whose inner measnre is discrete 
as above are stndied in Garr, Geman, Madan and Yor [2] with emphasis on hnancial 
application and called GGMY processes.) Observe that sample paths of fTSm look like 
their backgronnd driving Levy process as H is closer to l/a, while the dependence range 
gets longer and fTSm paths behave milder for greater H. 

For the reader’s convenience for comparison, we draw in Fignre El daily time series of 
TOYOTA shares on the Tokyo Stock Exchange, together with fTSm drawn in Figured 
It is observed that in short time the time series looks like fTSm with (0.8,1.6,P 2 ), while 
behaving in a Gaussian manner in longer time. 

To hnish this study, let us mention that the long time behavior result provides yet 
another way to simulate sample paths of fBm. For simplicity, consider a symmetric inner 
measure with a very simple structure, e.g. p{dx) = 2~^{S^i{dx) + Si{dx)), which reduces 
the random sequence to a sequence of iid Rademacher random variables 

Observe that 




i=l 


00 



Glearly, for sufficiently large h, the right hand side of the above behaves like 



i=l 
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{H, a, p) = (1.6, 0.7, pi) {H, a, p) = (1.0,1.2, ps) 



{H,a,p) = (0.8,1.6, P 2 ) 



{H,a,p) = (0.6,1.9, pi) 


Figure 1: Typical sample paths of fTSm (thick line) and of its background driving TS 
process (thin line) generated via the series representation. 




Figure 2: fTSm {H,a,p) = (0.8,1.6, P 2 ) (left thick) and (0.6,1.9, pi) (right thick) with 
(scaled) daily time series of TOYOTA shares on the Tokyo Stock Exchange; 247days (left 
thin) and 493days (right thin) up to 02/11/2005. 
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Theorem 16.41 (ii) tells us that this stochastic process (on [0,T]) approximates fBm. In 
order that its second moment is equal to that of a standard fBm, we set h = since 

then 


N 




2=1 


= h-^N 


2 a 

2 + a' 


_^2G ^ ^2G_ 


Therefore, for sufficiently large iV, the stochastic process 


2 a 

2 + a 


-N 


-1/2 N 


2=1 


iKH,a{i,Ti) : t G [0,T] 


can be used for simulation of standard fBm. 
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